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ABSTRACT 

Aims. We aim to constrain new starburst/AGN models of IRAS bright galaxies via their spectral energy distribution from the near- 
infrared to the radio. To this end, we determine the radio spectra for a sample of 31 luminous and ultraluminous IRAS galaxies 
(LIRGs/ULIRGs). 

Methods. We present here new high frequency VLA observations at 22.5 GHz and 8.4 GHz and also derive fluxes at other radio 
frequencies from archival data. Together with radio data from the literature, we construct the radio spectrum for each source. In 
the selection of data we have made every effort to ensure that these fluxes neither include contributions from nearby objects, nor 
underestimate the flux due to high interferometer resolution. 

Results. From our sample of well-determined radio spectra we find that very few have a straight power-law slope. Although some 
sources show a flattening of the radio spectral slope at high frequencies, the average spectrum shows a steepening of the radio spectrum 
from 1.4 to 22.5 GHz. This is unexpected, because in sources with high rates of star formation, we expect flat spectrum, free-free 
emission to make a significant contribution to the radio flux at higher radio frequencies. Despite this trend, the radio spectral indices 
between 8.4 and 22.5 GHz are flatter for sources with higher values of the far-infrared (FIR)-radio flux density ratio, q, when this is 
calculated at 8.4 GHz. Therefore, sources that are deficient in radio emission relative to FIR emission (presumably younger sources) 
have a larger thermal component to their radio emission. However, we find no correlation between the radio spectral index between 
1.4 and 4.8 GHz and q at 8.4 GHz. Because the low frequency spectral index is affected by free-free absorption, and this is a function 
of source size for a given mass of ionized gas, this is evidence that the ionized gas in ULIRGs shows a range of densities. 
Conclusions. The youngest LIRGs and ULIRGs are characterized by flatter average radio spectral indices from 1.4 to 22.5 GHz, 
and by a larger contribution to their high frequency, radio spectra from free-free emission. However, the youngest sources are not 
those that have the greatest free-free absorption at low radio frequencies. The sources in which the effects of free-free absorption are 
strongest are instead the most compact sources. Although these have the warmest FIR colours, they are not necessarily the youngest 
sources. 
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1. Introduction 

The nature of the power source for ULIRGs remains a much de- 
bated issue. The presence of either an AGN or a starburst in a 
given source is not evidence that one or the other is the principal 
power source for the infrared luminosity. The observation that 
ULIRGs fall on the same FIR-radio correlation as star-forming 
galaxies (e.g. Sopp, Alexander & Riley, 1990) was taken as ev- 
idence that these sources are powered predominantly by star 
formation. Analogous correlations exist between the Bry lumi- 
nosity and far-infrared luminosity (Goldader et al. 1997) and 
dense molecular gas mass and far-infrared luminosity (Gao & 
Solomon, 2004). However, evidence that AGN activity is not in- 
dependent of starburst activity have made such conclusions less 
secure. Farrah et al. (2003) show, in fact, that AGN and starburst 
luminosities are correlated over a wide range of IR luminosities. 
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As ULIRGs may commonly host both an AGN and a starburst, 
their relative contributions need to be quantified. 

Apart from the rarity of ULIRGs (there is no example closer 
than Arp 220 at 72 Mpc) their study is made difficult by the 
large and uncertain extinctions toward their centres. In many 
cases the extinction is not even negligible in the mid-infrared; 
it is likely that Ha emission detected from a ULIRG does not 
originate from the deeply embedded regions where the FIR lu- 
minosity is generated. The less extinguished, outer regions of a 
ULIRG may host a certain rate of star formation that causes a 
certain amount of Ha emission but this could be quite indepen- 
dent of whether it hosts a central starburst or an AGN. Murphy 
et al. (2001) find very few broad near-infrared lines in ULIRGs 
that would be the signature of an AGN, but note that very high 
extinctions could exclude their detection even at 2 /mi. Flores et 
al. (2004) recently showed that even when corrected by Hy data, 
the Ha luminosity can underestimate the star formation rate by a 
factor of 2 for ULIRGs. Both Goldader et al. (1995) and Valdes 
et al. (2005) find that near-infrared recombination lines (Paa 
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and Bry) are under-luminous in ULIRGs compared with what 
would be expected from a starburst of similar bolometric lumi- 
nosity. Their results indicate that even near-infrared observations 
may not penetrate the most obscured regions in these sources. 
This picture is consistent with the conclusions of Poggianti et 
al. (2001) that the extinction in the centre of ULIRGs is a func- 
tion of the age of the star formation episode, with younger stars 
being more heavily extinguished, as introduced by Silva et al. 
(1998) to describe star forming regions. The lack of correlation 
between super star clusters and HII regions found in luminous 
infrared galaxies by Alonso-Herrero et al. (2002) also supports 
a scenario of age-dependent extinction. In order to estimate the 
star formation rate in these sources tracers that do not suffer ex- 
tinction are preferred. 

Perhaps the most direct way to determine the power source 
in ULIRGs is via radio recombination lines that trace the ion- 
ized gas without suffering extinction. Observations of several 
lines between 1 .4 and 207 GHz of Arp 220 have been used by 
Anantharamaiah et al. (2000) to show that this source is powered 
by star formation. Unfortunately, very few ULIRGs are accessi- 
ble observationally to these kind of observations with existing 
instruments. 

One approach which has been used to distinguish between 
AGN and starburst power sources in ULIRGs has been to search 
for very compact radio continuum emission towards the nucleus 
or nuclei. Interestingly, radio sources of similar physical sizes 
have been identified as both AGN and compact starbursts, de- 
pending on the resolution of the observations. Nagar et al. (2003) 
use 15 GHz radio continuum data with a resolution of 150 mas 
to investigate the nature of 83 ULIRGs. One argument used by 
these authors to conclude that most ULIRGs are AGN powered 
is the compactness of the radio sources detected. The resolution 
of their data corresponds to 420 pc at the median redshift of their 
sample. However, the supernovae (SNe) detected in the NW nu- 
cleus of Arp 220 (Smith et al., 1998; Lonsdale et al. 2006) are 
within a region 0'.'2 x 0'.'4 (75 x 150 pc). Both of these studies 
found that no AGN is necessary to explain the IR luminosity of 
Arp 220. As long as the brightness temperature does not exceed 
10 6 K the compactness of radio nuclei alone does not support an 
AGN hypothesis. 

Rather than look at the radio morphology, we examine the 
spectrum of the radio emission. In star-forming galaxies the ra- 
dio spectrum is made-up of two components: a non-thermal, 
synchrotron component and a thermal, 'Bremsstrahlung' com- 
ponent, often referred to as 'free-free' emission. At GHz fre- 
quencies, the synchrotron spectrum, S v oc y ff , has a typical 
power-law slope, a ~ -0.8, while the free-free emission has 
a much flatter slope, a = -0.1. Free-free emission thus makes 
a larger contribution at higher radio frequencies, so that the ra- 
dio spectrum flattens at high frequencies. Half of the radio flux 
is expected to be of thermal origin at ~ 20 GHz (Condon, 
1992). Towards lower radio frequencies (below 1 GHz in non- 
starburst galaxies) the free-free optical depth becomes large and 
the (synchrotron-dominated) radio spectrum shows a sharp de- 
cline due to free-free absorption (see Fig. 4 of Condon, 1992). 
In sources with very intense star-formation, free-free absorption 
may be expected to occur at higher frequencies than in quies- 
cent objects; the 1.4 GHz fluxes from such objects are probably 
affected. 

In principle, a measure of the thermal radio flux from a 
ULIRG would be an excellent measure of the star formation 
rate, because, being emitted by the same gas from which re- 
combination lines originate, it traces gas ionized by young stars. 
Absorption of ionizing photons by dust may complicate this pic- 



ture slightly (Valdes et al, 2005) but the main difficulty is in esti- 
mating the thermal fraction. Although the fractional contribution 
from thermal emission to the total radio flux increases with fre- 
quency, the radio emission from ULIRGs is dominated by syn- 
chrotron radiation even at frequencies above 15 GHz. As we will 
see later, the radio spectral indices around 15 GHz are normally 
much steeper than the -0. 1 expected from a purely thermal spec- 
trum. The motivation for the 22 GHz observations described here 
was to better constrain the thermal radio flux. 

Because the thermal radio flux provides an estimate of 
the present star-formation rate, it is important in fitting star- 
burst/AGN models to the spectral energy distribution (SED) 
from the near-infrared to the radio. Bressan, Silva & Granato 
(2002) showed that deviations from the FIR/radio correlation 
could be used to derive the evolutionary status of a compact star- 
burst. Such deviations are expected in bursts of short duration be- 
cause at early times (a few xlO 6 yr) even the most massive stars 
formed in the burst will not have ended their lives as supernovae 
and no excess synchrotron emission should result from the star- 
burst activity. Such young sources then, should have a FIR/radio 
ratio above the mean FIR-radio correlation and should have flat- 
ter radio slopes due to the greater ratio of thermal/synchrotron 
emission. In Prouton et al. (2004) we illustrated the power of this 
technique with special emphasis on high frequency radio data. A 
second paper, Vega et al. (2007), hereafter paper 2, will be con- 
cerned with the model fits to the SEDs from the near-infrared to 
the radio. 

Here we present new radio data at 22.5 and 8.4 GHz, as 
well as archival radio data which has been re-reduced, and radio 
fluxes from the literature. The result is a set of 31, well-sampled 
radio spectra for infrared-luminous galaxies. In §[2] we describe 
the sample selection and in § [3] the data acquisition and reduc- 
tion. In § 2] we describe the radio spectra in detail, and in § [5] 
the implications of these spectra when compared to far-infrared 
fluxes and the FIR-radio flux ratio, q. Conclusions are reached 
in § [6] In an appendix we comment on individual objects and 
present the radio spectra in tabular form. 

2. Sample selection 

We have based our sample on that of Condon et al. (1991a) 
who made 8.4 GHz observations of the 40 ULIRGs brighter than 
5.25 Jy at 60 /mi in the IRAS Bright Galaxy Sample. In Prouton 
et al. (2004) we observed 7 of the ULIRGs observed by Condon 
et al. but restricted our attention to objects with a ratio of FIR to 
1 .4 GHz luminosity q > 2.5, in an attempt to preferentially select 
young starbursts. Here we have relaxed this restriction on q so as 
to maximize the sample size and include more diverse objects, 
and have observed a further 18 sources at 22.5 GHz. The re- 
maining sources were either too weak for observation at 22 GHz 
or were too confused with nearby sources to construct a reliable 
SED. We also re-observed 7 sources at 8.4 GHz for which the 
high resolution observations of Condon et al. (1991a) could not 
provide reliable integrated fluxes. Our final sample is therefore 
not a far-infrared flux limited sample like that of Condon et al., 
but simply contains all those brighter than 5.25 Jy at 60 fim for 
which an integrated 22.5 GHz flux could be reliably obtained. 

In addition to the new high frequency observations, we have 
obtained, either from the literature or from reduction of archival 
VLA data, radio fluxes at as many frequencies as possible. In 
total, of the 40 ULIRGs of the sample of Condon et al., 31 have 
radio data at 3 or more frequencies and 29 of these have fluxes 
at 22.5 GHz. 

All our sample galaxies are at a redshift z < 0.1. 
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3. Observations and data reduction 

3.1. New 22.5 and 8.4 GHz data 

22.5 GHz observations were carried out on 2004, August 13 and 
14 with the VLA in D-configuration. Total integration times for 
the sources were between 16 and 24 minutes, with scans on- 
source being interleaved between scans of phase calibrators ev- 
ery 8 minutes. Before each new source was observed, a pointing 
scan was carried out on the phase calibrator at 8.4 GHz in order 
to determine a pointing correction. 3C 286 was observed as an 
absolute flux calibrator on 13 August and 3C 48 on 15 August. 

The weather conditions on August 13 were not ideal for high 
frequency radio observations, data were taken in the presence 
of thunder storms. As water vapour affects the measured flux 
for a given source at 22 GHz, and the flux calibrator may have 
been observed through a different cloud density than some of 
the phase calibrators, the bootstrapping of the phase calibra- 
tor fluxes to those of the flux calibrators is not guaranteed to 
work. After flux calibration based on the clean component mod- 
els for 3C 286 and 3C 48 provided by NRAO the phase calibra- 
tor fluxes were within 20% of those given by the VLA calibrator 
database. Only in the case of 1310+323 (the phase calibrator for 
UGC 8387) was there a larger discrepancy (a factor of 2 less than 
given in the database). 

Maps were made using 'natural' weighting in the uv -plane 
and in the case of UGC 6436, IZW 107 and IRAS 14348-144 a 
uv-taper of 30 kA was applied to avoid the possibility of missing 
flux due to extended but weak structure. The final maps, shown 
in Fig.Q] had a typical resolution of 4" (6" where tapered) and 
after CLEANing had a typical rms noise level of 0.12 mJy/beam. 
Fluxes were measured by direct integration on the maps. These 
fluxes (reported in table [1} have errors which include the rms 
noise and the flux calibration errors described above. 

8.44 GHz observations were made on 2005, July 12, 20 and 
27 with the C-configuration of the VLA. Integration times of 1 1 
to 14 minutes were used and interleaved with observations of 
nearby phase calibrators. 3C 286 and 3C 48 were used as abso- 
lute flux calibrators and, as for the 22 GHz data, calibration made 
use of clean component models for these sources. Derived fluxes 
for all phase calibrators agreed with those in the VLA calibrator 
database to better than 5%. All sources were also self-calibrated 
to correct for phase only. Maps were made by natural weighting 
the uv-data and CLEANed in the standard way. Resolution and rms 
in the resulting maps were typically 3" and 0.03 mJybeam -1 . 
The radio maps for these sources are presented in Fig.|2]and the 
corresponding integrated fluxes in table [2] 

3.2. Archival VLA data 

In addition to the above data and fluxes taken from the literature 
the VLA archive was searched for data at additional frequencies 
so as to maximize the coverage of the radio spectrum. Data suit- 
able for the measurement of integrated fluxesj were downloaded 
and reduced following standard procedures. Where the sources 
were of sufficient strength, a single iteration of self-calibration 
was applied allowing for the correction of antenna phases only. 
Fluxes were measured by direct integration on the CLEANed maps 
and are presented in tabled For UGC 08058 (Mrk 231) archival 
data were downloaded and reduced, despite the presence of mea- 
sured fluxes in the literature, because of evidence that the radio 



i.e., where the source extent was smaller than the largest angular 
scale which can be mapped by the VLA in the given configuration and 
frequency band. 



Table 2. 8.4 GHz fluxes for the newly observed objects. 



Source 


Source flux 


Phase cal 


Phase cal flux 




(mJy) 


(J2000) 


(Jy) 


IC 1623 


54.6 ± 3.0 


0118-216 


0.854 ± 0.002 


Arp 148 


11.8 + 0.6 


1130+382 


1.338 + 0.005 


UGC 6436 


5.2 ±0.3 


1122+180 


0.639 ± 0.003 


IRAS 14348-1447 


10.8 + 0.5 


1439-155 


0.509 ± 0.003 


IZw 107 


14.3 ± 0.7 


1500+478 


0.341 ±0.001 


NGC 6286 


31.5 ± 1.6 


1638+573 


1.261 +0.003 


Mrk 331 


22.3 ± 1.0 


0010+174 


0.670 ± 0.002 



emission from this source is variable. A significant fraction of 
its radio flux must therefore come from an AGN. The 5 fluxes 
listed in table [3]were all measured on the same date. The maps 
of sources reduced from the VLA archive are shown in Fig. [3] 

All calibration and mapping was carried using the NRAO 
Astonomical Image Processing System (aips). The radio fluxes 
for our sample galaxies (in the range 100 MHz to 100 GHz) are 
listed in table 1 in the appendix, along with notes on selected 
sources. 

3.3. Data at other frequencies 

In order to make comparisons between our radio data and in- 
frared data, we also collected J,H and K-band fluxes from the 
2MASS Extended Source Catalogue (Jarrett et al., 2000) and 12, 
25, 60 and 100 fim fluxes from the IRAS Faint Source Catalogue 
(Moshiretal., 1990). 

In the collection of fluxes at other frequencies from the liter- 
ature, care was taken to ensure that the fluxes were not confused 
or, in the case of interferometric radio data, underestimated due 
to high resolution. While most sources are compact and do not 
have other strong sources nearby, a small number of sources re- 
quired corrections to some fluxes. 

Confusion is most likely to be a problem for the IRAS fluxes 
due to the very low resolution. To try to quantify any possible 
problems with confusion in the absence of a higher resolution 
FIR survey, we used the NVSS and the far infrared-radio cor- 
relation. In those cases where an unrelated radio source was 
within the IRAS beam we reduced the IRAS flux by the frac- 
tional contribution of this source in the radio. Few sources were 
affected, but one source (IRAS 03359+1523) has its FIR emis- 
sion confused with a source of similar strength. Because the 
error on a correction of this scale would render the fluxes un- 
reliable this source will not be considered in the model fitting 
described in paper 2. The ULIRGs for which the IRAS fluxes 
were reduced because of the presence of a nearby contaminating 
source were as follows: UGC 6436 (26%), IRAS 12112+0305 
(15%), NGC 6286 (6%), UGC 5101 (8%). 

The full list of fluxes for the sample galaxies from the near- 
infrared to the radio will be presented in paper 2, where they are 
used to model the spectral energy distributions. 
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Table 1. 22.5 GHz fluxes derived from the new VLA observations. * Observations described in Prouton et al. (2004) 



Source 


R.A. 


Dec. 


Source flux 


Phase cal 


Phase cal flux 




(J2000) 


(J2000) 


(mJy) 


(J2000) 


(Jy) 


NGC 34* 


00 11 06.550 


-12 06 26.33 


7.41 +0.14 


2358-103 


0.5719 ±0.005 


IC 1623 


01 07 47.49 


-17 30 26.3 


21.7 + 2.0 


0050-094 


0.805 


+ 0.007 


CGCG436-30* 


01 20 02.722 


+ 14 21 42.94 


8.73 ±0.31 


0121+118 


1.177 


± 0.007 


IRAS 01364-1042* 


01 38 52.870 


-10 27 11.70 


3.97 ±0.19 


0141-094 


0.7395 


± 0.003 


IIIZw 35 


01 44 30.53 


+ 17 06 08.5 


9.7 ±0.3 


0121+118 


1.449 


±0.011 


UGC 2369 


02 54 01.81 


+ 14 58 14.5 


5.3 ±0.3 


0238+166 


2.041 


±0.015 


IRAS 03359+1523 


03 38 47.12 


+ 15 32 53.5 


4.5 ±0.3 


0325+226 


0.768 


± 0.006 


NGC 1614 


04 34 00.02 


-08 34 44.98 


21.0 ± 1 


0423-013 


7.480 


± 0.058 


NGC 2623 


08 38 24.08 


+25 45 16.40 


18.3 ±0.2 


0830+241 


1.822 


±0.013 


IRAS 08572+3915* 


09 00 25.390 


+39 03 54.40 


3.18 ±0.3 


0927+390 


7.688 


±0.15 


UGC 4881 


09 15 55.52 


+44 19 57.5 


3.5 ±0.5 


0920+446 


2.140 


±0.010 


UGC 5101 


09 35 51.72 


+61 21 11.29 


18.0 ± 1.0 


0958+655 


0.972 


±0.010 


IRAS 10173+0828 


10 20 00.22 


+08 13 34.5 


2.9 ±0.4 


1041+061 


1.672 


± 0.009 


Arp 148 


11 03 53.95 


+40 50 59.5 


6.1 ±0.3 


1146+399 


1.330 


±0.012 


UGC 6436 


11 25 45.15 


+ 14 40 29.4 


2.8 ±0.5 


1118+125 


1.238 


±0.016 


IRAS 12112+0305 


12 13 46.1 


+02 48 41.5 


5.7 ±0.3 


1224+035 


0.788 


± 0.006 


UGC 8387 


13 20 35.34 


+34 08 22.19 


7.5 ± 2.0 


1310+323 


0.963 


± 0.042 


UGC 8696 


13 44 42.11 


+55 53 13.15 


14.3 ±0.3 


1419+543 


0.935 


± 0.023 


IRAS 14348-1447 


14 37 38.30 


-15 00 24.0 


4.5 ±0.5 


1448-163 


0.376 


± 0.006 


IZw 107 


15 18 06.07 


+42 44 45.19 


3.6 ± 0.3 


1521+436 


0.195 


± 0.005 


IRAS 15250+3609 


15 26 59.40 


+35 58 37.53 


4.2 ±0.3 


1521+436 


0.195 


± 0.005 


NGC 6286 


16 58 31.70 


+58 56 14.50 


12.2 ± 1.5 


1638+573 


1.259 


±0.018 


NGC 7469* 


23 03 15.623 


+08 52 26.39 


17.50 ±0.50 


2320+052 


0.5661 ± 0.004 


IC 5298* 


23 16 00.690 


+25 33 24.20 


3.86 ±0.18 


2321+275 


0.7178 ±0.009 


Mrk331* 


23 51 26.802 


+20 35 09.87 


9.95 ±0.29 


2358+199 


0.2683 


± 0.002 



Table 3. Fluxes derived from archival VLA data. * No flux calibrator observed. Flux of 0513-219 set to this value using the VLA 
calibrator database. The last column, n, shows the scaling of the contours in Fig. [3] 



Source 


Frequency 


Flux 


Phase cal 


Phase cal flux 


Program ID 


Array 


Date 


n 




(GHz) 


(mJy) 


(J2000) 


(Jy) 










NGC 34 


4.8 


28.4 ± 0.3 


2358-103 


0.8080 ±0.0013 


AM451 


B 


1994-M-24 


6 


IIIZw 35 


15 


10.0 ± 1 


0149+189 


0.3538 ± 0.003 


TYP04 


B 


2002-sep-06 


-1 


NGC 1614 


8.4 


40.7 ± 0.3 


0423-013 


3.1178 ±0.004 


AK331 


C 


1993-aug-03 


6 


IRAS 05189-2524 


15 


7.8 ±0.2 


0513-219 


0.92 ± 0.01* 


AC624 


C 


2003-jan-05 


4 


IRAS 10173+0828 


4.8 


5.88 ± 0.2 


1150-003 


1.7056 ±0.004 


AK184 


CD 


1988-mar-28 


4 




15 


3.1 ±0.3 


1150-003 


0.9135 + 0.004 


AK184 


CD 


1988-mar-28 


1 


Arp 148 


4.8 


14.0 ±0.5 


1130+382 


0.7813 ±0.002 


AH333 


A 


1989-jan-09 


5 




15 


8.0 ± 1 


1130+382 


0.6237 + 0.007 


AJ105 


C 


1984-may-24 


4 


Arp 299 


8.4 


243 ±8 


1219+484 


0.7873 + 0.003 


AW641 


C 


2005-aug-20 


-1 




15 


160 ±3 


1219+484 


0.8489 ± 0.012 


AS568 


c 


1999-jan-21 


2 




22.5 


74 ±2 


1128+594 


0.4509 ± 0.004 


AN095 


c 


2001-aug-25 


2 


IRAS 12112+0305 


4.8 


14.7 ±0.2 


1150-003 


1.7056 ±0.004 


AK184 


CD 


1988-mar-28 


3 




15 


5.7 ±0.2 


1150-003 


0.9135 + 0.004 


AK184 


CD 


1988-mar-28 


1 


UGC 08058 


1.4 


274 ±3 


1219+484 


0.5987 + 0.001 


BU006Y 


B 


1995-nov-17 


1 




4.8 


265 ± 1 


1219+484 


0.7885 ± 0.002 


BU006Y 


B 


1995-nov-17 


3 




8.4 


189 ± 1 


1219+484 


0.9363 ± 0.002 


BU006Y 


B 


1995-nov-17 


4 




15 


146 ± 1 


1219+484 


1.0916 ±0.005 


BU006Y 


B 


1995-nov-17 


3 




22.5 


137 ±2 


1219+484 


1.1123 ±0.008 


BU006Y 


B 


1995-nov-17 





UGC 8387 


15 


12 ±2 


3C286 


3.45375 


AN35 


C 


1985-aug-30 





NGC 5256 


8.4 


24.1 ± 1 


1417+461 


0.5880 ± 0.003 


AP233 


C 


1992-apr-28 


4 




15 


21.0 ± 1 


1349+536 


0.6003 ± 0.006 


AM 198 


CD 


1987-jan-22 


2 


UGC 8696 


15 


29.5 ± 1 


1349+536 


0.6003 ± 0.006 


AM198 


CD 


1987-jan-22 


2 


NGC 7469 


8.4 


50.2 ± 0.6 


2320+052 


0.5897 ± 0.05 


AP233 


C 


1992-apr-28 


3 



Table 4. Far-infrared-radio flux ratios, q, at 1.4 and 8.4 GHz and radio spectral indices for the whole sample in various frequency 
ranges. The value of FIR = 1.26 x 1(T 14 (2.58 S(,q + S ioo), as used to calculate q, is also given. 



■Tame 


rlK 

(/10 W m ) 




98.4 


/T 1 ' 4 

a 22.5 


^,1.4 
a 4.8 


a 84 


^,4.!! 


„,8.4 




n/8.4 
"22.5 


A- 15 

u 22.5 


JGC 34 


7.69 ± 0.02 


/-> Af\ , A AO 

2.49 ± 0.02 


O 1 O i A AO 

3.13 ± 0.02 


a an i a ao 

-0.79 ± 0.02 


a an i a ao 

-0.70 ± 0.03 


A O O i A A A 

-0.83 ± 0.04 


1 1 O i A 1 A 

-1.12 ± 0.10 






A TO 1 A f\£L 

-0.73 ± 0.06 




C 1623 


11.2 ± 0.63 


2.08 ± 0.03 


2.74 ± 0.03 


A O O . A A A 

-0.88 ± 0.04 


A A 1 1 

-0.77 ±0.11 


A OC A AO 

-0.85 ± 0.02 


1 A 1 A O O 

-l.Ol ± 0.23 


~ 




A A A A 1 A 

-0.94 ± 0.10 


~ 


,GCG 436-30 


A l ~If\ A AO 

4.70 ± 0.03 


r\ ,i f\ A A1 

2.40 ± 0.01 


O AA . A AO 

2.99 ± 0.02 


A f O . A AO 

-0.63 ± 0.02 


A /A . A AC 

-0.69 ± 0.05 


A T1 i A AO 

-0.77 ± 0.03 


A A A . A 1 O 

-0.94 ± 0.12 






A O O 1 A C\£L 

-0.38 ± 0.06 




RAS 0136-1042 


O AO A AO 

3.02 ± 0.02 


2.71 ± 0.02 


O AA A AO 

2.99 ± 0.02 


A CA A AO 

-0.50 ± 0.02 


A O A .A AC 

-0.24 ± 0.05 


A O a . A A A 

-0.37 ± 0.04 


-0.65 ± 0.12 






a a a .a Aa 

-0.74 ± 0.07 




IIZw 35 


c i—i a a or 

5.74 ± 0.35 


r\ r O A AO 

2.58 ± 0.03 


O OA . A AO 

2.89 ± 0.03 


A CO . A AO 

-0.52 ± 0.02 


A O C A A C 

-0.35 ± 0.05 


f\ \ A . A AO 

-0.40 ± 0.03 


A CO A 1 O 

-0.52 ± 0.13 


-1.17 ± 


0.19 


A ao A A/" 

-0.72 ± 0.06 


A AO AO/' 

-0.08 ± 0.26 


JGC 2369 


3.94 ± 0.16 


2.33 ± 0.02 


O AA A AO 

2.90 ± 0.03 


A OA A AO 

-0.80 ± 0.02 




a ao A AO 

-0.73 ± 0.03 








A AO A AO 

-0.93 ± 0.08 




RAS 03359+1523 


r\ o O . A 1 A 

2.88 ± 0.14 


r\ s f\ r\ z\o 

2.60 ± 0.03 


O O A . A AO 

2.84 ± 0.03 


A CO A AO 

-0.53 ± 0.03 


A A 1 .A AC 

-0.41 ± 0.05 


A OO . A AO 

-0.32 ± 0.03 


-0.11 ± 0.13 






A A 1 .A AO 

-0.91 ± 0.08 




JGC 1614 


-\ A £L r\ ac 

14.6 ± 0.45 


o AC f\ rvO 

2.45 ± 0.02 


O AO i A A 1 

2.98 ± 0.01 


a /n , A AO 

-0.68 ± 0.02 


A £. O i A 1 A 

-0.63 ± 0.14 


A /o , A AO 

-0.68 ± 0.02 


A a O i A O 1 

-0.78 ± 0.31 






a z:a 1 a ac 

-0.67 ± 0.05 




pj i n nr 1 OH OCO /I 

RAS 05189-2524 


C OA i A in 

5.89 ± 0.19 


O l-l A , f\ AO 

2.74 ± 0.02 


O 1 A ■ A AO 

3.14 ± 0.03 




A A O ■ A AC 

-0.43 ± 0.05 


A CO i A AO 

-0.52 ± 0.03 


A a 1 i A 1 o 

-0.71 ± 0.13 


-0.65 ± 


0.10 






JGC 2623 


11.0 ± 0.35 


A ^A , f\ /"VO 

2.49 ± 0.02 


O AO i A AO 

2.92 ± 0.03 


A /A . A A1 

-0.60 ± 0.01 


A O A i A 1 A 

-0.39 ± 0.14 


A ff . f\ AO 

-0.55 ± 0.03 


A A 1 ■ A O O 

-0.91 ± 0.32 






A ti a , A AC 

-0.67 ± 0.05 




n ao n O C a O . T A 1 c 

RAS 08572+3915 


r\ 1-1—7 j-\ AO 

2.97 ± 0.02 


O Oa A A/1 

3.27 ± 0.04 


O OA A AO 

3.29 ± 0.02 


A 1 1 .A AC 

-0.11 ± 0.05 


A AO i A AO 

0.02 ± 0.08 


A AO . A A/ 

-0.03 ± 0.06 


A 1 O i A 1 O 

-0.13 ± 0.12 






A Oj' A 1 1 

-0.26 ±0.11 




JGC 4881 


3.22 ± 0.13 


2.37 ± 0.02 


O AA A AO 

2.99 ± 0.03 


A O C A A C 

-0.85 ± 0.05 




A OA . A AO 

-0.80 ± 0.03 








f\ f\ A A 1 C 

-0.94 ± 0.15 




JGC 5101 


Z" O O A 

6.28 ± 0.29 


1 AA f\ AO 

1.99 ± 0.03 


O CA A AO 

2.50 ± 0.03 


A O 1 .A AO 

-0.81 ± 0.02 


A f A A 1 'I 

-0.69 ± 0.17 


A f f . A AO 

-0.66 ± 0.03 


A CA . A O O 

-0.59 ± 0.38 


-0.97 ± 


0.58 


1 AA . A AO 

-1.09 ± 0.08 


1 O f A O O 

-1.26 ± 0.83 


RAS 10173+0828 


2.58 ± 0.14 


O O -1 f\ AC 

2.84 ± 0.05 


O 1 A . A AO 

3.10 ± 0.03 


A A A .A A/ 

-0.44 ± 0.06 


A A O A AO 

-0.42 ± 0.08 


A O A .A A/" 

-0.34 ± 0.06 


-0.15 ± 0.11 


-0.96 ± 


0.19 


-0.63 ± 0.15 


A 1 f A A O 

-0. 16 ± 0.42 


RAS 10565+2448 


c a a r\ oo 

5.84 ± 0.22 


f\ A A , r\ AO 

2.44 ± 0.02 


O A A . A AO 

3.04 ± 0.03 




A T~t . A AO 

-0.77 ± 0.03 


A ao . A AO 

-0.78 ± 0.03 


A O 1 .A AA 

-0.81 ± 0.09 










s.rp 148 


O O O A 1 o 

3.23 ± 0.13 


o oa a ao 

2.37 ± 0.02 


O OT A A A 

2.87 ± 0.04 


r\ /~ a , r\ r\o 

-0.64 ± 0.02 


-0.78 ± 0.04 


A f O . A AC 

-0.63 ± 0.05 


-0.31 ± 0.16 


-0.66 ± 


0.26 


A /'a A 1 A 

-0.67 ± 0.10 


A z' a A 

-0.67 ± 0.33 


JGC 6436 


o o o .a nz' 

3.23 ± 0.06 


-\ /- r- z\ z\o 

2.65 ± 0.02 


O OO A AO 

3.22 ± 0.03 


A TA A A^7 

-0.70 ± 0.07 




A ao . A A A 

-0.73 ± 0.04 








A z' O . A 1 A 

-0.63 ± 0.19 




Vrp 299 


47.2 ± 2.13 


A/ A Al 

2.26 ± 0.03 


1 n I .A AO 

2.71 ± 0.02 


A OA j A AO 

-0.80 ± 0.02 


A A O A 1 A 

-0.43 ± 0.10 


A CO j A AO 

-0.58 ± 0.03 


A AA . A O 1 

-0.90 ± 0.21 


-0.72 ± 


0.07 


1 1 A A A 

-1.21 ± 0.04 


1 AA A AO 

-1.90 ± 0.08 


n a o irti 10 i none 

RAS 12112+0305 


4.02 ±0.17 


O ZT ZT , f\ 

2.66 ± 0.02 


O AO ■ A AO 

3.03 ± 0.03 


A C 1 i A AO 

-0.51 ± 0.02 


A on , f\ AO 

-0.37 ± 0.03 


f\ A a , A AO 

-0.47 ± 0.03 


r\ /~{~\ i A A A 

-0.69 ± 0.09 


-0.97 ± 


0.11 


a ca a Aa 

-0.57 ± 0.07 


A AA , A 1 / 

-0.00 ± 0.16 


Ty— i oaco 

JGC 8058 


14.2 ± 0.54 


O 1 /I .A A^ 

2.14 ± 0.02 


r\ OA A AO 

2.30 ± 0.02 


A o c f\ A 1 

-0.25 ± 0.01 


A AO A A 1 

-0.03 ± 0.01 


A O 1 A A 1 

-0.21 ± 0.01 


A S A A A 1 

-0.60 ± 0.01 


-0.45 ± 


0.01 


A O O A AO 

-0.33 ± 0.02 


-0.16 ± 0.04 


JGC 8387 


O 1 O i A A A 

8. 18 ± 0.49 


T , A AO 

2.32 ± 0.03 


O OA ■ A AO 

2.80 ± 0.03 


A AC ■ A 1 A 

-0.95 ± 0.10 


a ao i A AC 

-0.78 ± 0.05 


A £L 1 ■ A AO 

-0.61 ± 0.03 


A O A i A 1 O 

-0.24 ± 0.13 


-1.84 ± 


0.30 


1 C£L 1 A O O 

-1.56 ± 0.28 


-1.16 ± 0.78 


JGC 5256 


3.78 ± 0.17 


1 AA A AO 

1.90 ± 0.03 


O s r\ A AO 

2.62 ± 0.03 


-0.77 ± 0.04 


a a a a a /* 

-0.77 ± 0.06 


A AO A AO 

-0.93 ± 0.03 


1 O O A 1 A 

-1.28 ± 0.14 


-0.23 ± 


0.11 


c\ a a A 1 1 

-0.48 ±0.11 


a oc a oa 

-0.85 ± 0.27 


JGC 8696 


r\ nr a aa 

9.75 ± 0.30 


2.25 ± 0.02 


O TO . A AO 

2.78 ± 0.03 


A O O . A A 1 

-0.83 ± 0.01 


A C O A 1 a 

-0.58 ± 0.17 


A /~ a A AO 

-0.67 ± 0.03 


A O O . A O A 

-0.88 ± 0.39 


-0.67 ± 


0.10 


— 1.13 ± 0.06 


1 a A A 1 A 

-1.79 ± 0.10 




j. iz ± u. l^t 


A.j I ± U.UZ 


Z.oy ± VJ.VJH 


— u. 10 it u.\Jh 




A f\l -t- A A^ 
— U.U / ± U.UJ 








A SQ -1- A 1 ^ 




Zw 107 


4.34 ± 0.13 


T OA i A AO 

2.36 ± 0.02 


2.91 ± 0.03 


A AC i A AO 

-0.95 ± 0.03 


—0.3 / ± 0.03 


A Ht\ i A A/1 

-0. /0 ± 0.04 


1 A A i A 1 1 

-1.44 ± 0.13 






1 A A 1 A 1 1 

- 1 .40 ± 0. 1 1 




RAS 15250+3609 


3.11 ±0.12 


2.76 ± 0.02 


2.90 ± 0.03 


-0.45 ± 0.03 


-0.09 ±0.19 


-0.18 ±0.04 


-0.38 ± 0.42 


-0.47 ± 


0.65 


-0.93 ± 0.09 


-1.59 ±0.94 


up 220 


47.9 ± 1.41 


2.59 ± 0.02 


2.94 ± 0.03 


-0.46 ± 0.03 


-0.36 ± 0.08 


-0.44 ± 0.03 


-0.63 ±0.19 


-0.51 ± 


0.09 


-0.50 ± 0.08 


-0.49 ±0.17 


JGC 6286 


5.08 ±0.13 


1.94 ±0.02 


2.63 ± 0.02 


-0.95 ± 0.05 


-0.88 ±0.13 


-0.90 ± 0.03 


-0.93 ± 0.28 






-1.05 ±0.15 




JGC 7469 


13.3 ±0.08 


2.29 ± 0.01 


2.85 ± 0.01 


-0.84 ± 0.01 


-0.77 ±0.10 


-0.71 ±0.02 


-0.59 ±0.21 






-1.07 ±0.03 




C 5298 


4.46 ± 0.02 


2.53 ± 0.02 


3.16 ±0.02 


-0.79 ± 0.02 


-0.74 ± 0.05 


-0.81 ± 0.04 


-0.94 ±0.12 






-0.76 ± 0.07 




/IRK 331 


8.71 ±0.03 


2.52 ± 0.01 


3.02 ± 0.02 


-0.70 ± 0.02 


-0.67 ± 0.05 


-0.64 ± 0.03 


-0.59 ±0.12 






-0.81 ±0.05 




/lean 




2.441 ± 0.004 


2.906 ± 0.005 


-0.671 ± 0.007 


-0.521 ±0.018 


-0.596 ± 0.007 


-0.698 ± 0.040 


-0.790 ± 0.079 


-0.814 ±0.020 


-0.842 ±0.136 
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4. Results 

4.1. Radio spectral indices 

Before describing the observed radio spectra, we recall the ex- 
pected, general form of the radio spectrum for a compact star- 
burst, as outlined in §Q] At frequencies, v ~ 10 GHz, a steep 
synchrotron spectrum dominates, with S v oc V, a ~ -0.8. At 
both lower and higher frequencies the spectrum may be expected 
to flatten due to free-free absorption (low frequencies) and free- 
free emission (high frequencies). 

In fig. |4] we show the radio spectra for all our sources and ta- 
ble 2]lists the radio spectral indices, a, for the sample objects in 
several frequency ranges. Also shown are the far-infrared-radio 
flux ratios, q = log(FIR/3.15 x 10 12 Hz)/S v , where S v is the 
radio flux in units of Wm _2 Hz _1 (Condon et al., 1991a). We 
will use the notation, q\ 4, to refer to this quantity calculated at 
1.4 GHz. 

The mean value of a between 1.4 and 4.8 GHz (all here- 
after) is -0.521. This steepens steadily towards higher frequen- 
cies, with al l = -0.698, a* 22 \ = -0.813, and a^ 5 = -0.842 (al- 
though only 12 objects have values in this last range). This trend 
can also be seen in the lower-right panel of fig. [4] At low frequen- 
cies, the flatter spectra are probably due to free-free absorption. 
However, the fact that the spectra steepen to frequencies as high 
as 22.5 GHz is unexpected, since these sources should have very 
significant free-free emitting components, which would tend to 
flatten the spectra at high frequencies. Some sources even show 
spectral indices, ofjc, that are steeper than what is expected 
from pure synchrotron emission. The general trend, and these 
sources in particular, are further discussed in § 15.31 

In Fig. [5] we show how the radio spectral index is related to 
q. There is a clear tendency for higher values of q\ 4 to be asso- 
ciated with radio spectral indices which are flatter, on average, 
from 1.4 to 22.5 GHz. There is also a clear correlation between 
q\ 4 and a l A, with a much weaker correlation between q\^ and 
at *. On the other hand, 4 is not correlated with either cr^g or 
a^. This trend would be observed if the 1.4 GHz radio fluxes 
were reduced by the effects of free-free absorption. In this case, 
both gi 4 and a\ ■! are affected. The 4.8 GHz flux should be much 
less affected than the 1 .4 GHz flux because the free-free optical 
depth depends on frequency as v -21 . The 8.4 GHz flux is free 
from the effects of absorption, and so no correlation is seen be- 
tween ^8.4 and ff^-g or a^. 

At high frequencies we see a correlation between q$, .4 and 
a 22 5 ■ T ne corre l a ti on seen in the lower right panel of fig[5]shows 
that the fractional contribution of free-free emission is larger in 
sources with higher values of q^A- Therefore, despite the fact 
that the average source shows a radio spectrum which steepens 
towards higher frequencies, flatter spectral indices tend to occur 
in sources with higher q$, .4. This is expected if the synchrotron 
emission is powered by supernova explosions because of the de- 
lay between the formation of massive stars and the occurrence 
of the first supernovae. 

The trend for higher values of q to be associated with flatter 
radio spectra over the whole range from 1.4 to 22.5 GHz, ab£ 5 , 
will be used as a constraint on the starburst age in the model 
fitting of paper II. 

4.2. Possible AGN contribution to the radio fluxes 

If there are ULIRGs in our sample that contain an AGN it is pos- 
sible that this AGN contributes to the radio flux of the source. If 
the power-law index of the radio emission is different from that 



of the star formation we may expect to see a correlation between 
the radio spectral index and the near-infrared colours, which are 
sensitive to the presence of an AGN. In fig. [6] we plot the ra- 
dio spectral index against the near-infrared colour index (J-K). 
There is no indication of a general correlation, which suggests 
that most sources do not contain a radio loud AGN. However, 
we do note that UGC 8058 (Mrk 231) which has the highest 
value of (J-K) also has a very flat radio spectral index. As men- 
tioned above, there is evidence that the AGN in this source is 
radio loud. 

4.3. Radio spectral indices and IRAS colour, (f 60 / f m ) 

In fig. 13 we plot the radio spectral indices against the IRAS 
colour (/6o//ioo)- a \i shows a correlation with IRAS colour (in 
good agreement with the results of Sopp & Alexander, 1992) 
but this correlation almost disappears if a\ 2 5 is considered in- 
stead. Because warmer IRAS colours are accompanied by larger 
FIR fluxes (Young et al. 1989, Soifer et al. 1989) warmer IRAS 
colours imply higher rates of star formation. Therefore sources 
with warmer colours should have larger masses of ionized gas. 
However, the weakness of the correlation between (/eo/Zioo) and 
&22 5 suggests that this effect is not the main driver of the stronger 
correlation seen with a] g. 

The compactness of the regions is probably what causes the 
correlation in the upper panel of fig. Q If warmer IRAS colours 
were associated with more compact sources, then greater free- 
free absorption would be seen in warmer sources and a correla- 
tion between (/60//100) and a l 4 ^ would result, even if the mass of 
ionized gas does not change with (fm/fioo)- Source geometry is 
considered in more detail in § [5] 

5. Discussion 

5.1. The radio spectral index and the FIR-radio flux ratio, q 

In fig.|4]there are several sources in which both a low and a high 
frequency spectral flattening are seen (e.g. IRAS 12112+0305). 
This flattening is due to free-free absorption at low frequencies 
and free-free emission at high frequencies. However, there are 
also sources which show evidence of free-free emission at high 
frequencies, with no sign of free-free absorption at low frequen- 
cies (e.g.CGCG436) and others that show signs of free-free ab- 
sorption, with no sign of free-free emission (e.g. NGC 2623). 
The simplest explanation for these contrasting radio spectra is 
probably variations in the compactness of the emitting regions. 

For a region of radius r and density Nh maintained 
in ionization equilibrium by a source of ionizing radiation 
of fixed luminosity, the total ionization/recombination rate, 
(4/3)nT 3 A^Q'B, where ag is the recombination coefficient, is 
constant (Osterbrock, 1989). The free-free optical depth, Tff oc 
« 2 r, for a density of ionized gas n e . Therefore, as long as Tg 
remains below unity (certainly true at 22 GHz) the integrated 
free-free emission, oc Tff r 2 , from a region with fixed ionizing 
luminosity, is independent of the source size. 

However, where this gas absorbs a luminosity source, the 
absorption is proportional to e~ T " (foreground screen) and the 
dependence of Tff on source size means that the compactness 
of a region very strongly affects the amount of absorption. For 
a given ionizing luminosity and synchrotron luminosity, an ex- 
tended region may show no evidence of free-free absorption at 
1 .4 GHz, whereas a more compact region will. 

This effect is probably behind the the diverse spectral shapes 
of fig. [4] and is also important in understanding fig. [5] 
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Fig. 1. 22.5 GHz images of the newly observed ULIRGs in our sample. Contour levels begin at +0.5 mJy beam 1 and are spaced by 
a factor of 2'^ 2 . The synthesized beam is shown in the bottom left of each panel. 
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Fig. 2. Contour maps for the new 8.4 GHz observations. The synthesized beam is shown in the bottom left of each panel. NGC 6285 
is the companion of NGC 6286 and its radio map is included only for completeness. 



In considering the relation between the radio spectral indices emission is also proportional the star formation rate, but as the 

at various frequencies and q, shown in fig. [5] it is useful to keep emission is delayed by main sequence lifetime of massive stars, 

in mind how each are related to the star formation rate. Both it measures the star formation rate ~ 10 7 years prior to the cur- 

the FIR and the integrated free-free emission are expected to be rent epoch. Therefore, in objects characterized by bursts of star 

proportional to the current star formation rate. The synchrotron formation, such as those in our sample, we expect both a and q 
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Fig. 6. The relation between the near-infrared colour J-K and the Fig. 7. Top: Correlation between the radio spectral index be- 
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22.5 GHz (bottom). Objects at the extremes of the distributions Bottom: Same plot for the 8.4-22.5 GHz spectral index. Objects 
are labeled. at the extremes of the distributions are labeled. Note that that 

UGC 8058 may have a radio loud AGN. 



to be a function of the age of the starburst. Values of q are higher 
where the current star formation rate is greater than that ~ 10 7 
years ago, and lower where the star formation rate has declined. 
In this scenario, q is a measure of the "age" of a starburst, with 
higher values associated with younger objects. 

In fig. [5] as well as the correlations between q\n and a l 4 g, and 
between q%ji and we note the lack of a correlation between 

^8.4 and a l 4 g. If large values of q%^ are due to a delay in the 
production of synchrotron emission relative to that of thermal 
dust and free-free emission, as the correlation between q$ .4 and 
suggests, then the lack of correlation between gs.4 and ar^'g 
implies that the age of the starburst is not the only effect. 

Because source compactness strongly affects the amount of 
free-free absorption, but does not affect the free-free emission 
(for Tff <k 1) the presence of free-free emission at 22.5 GHz 
does not imply the presence of free-free absorption at 1 .4 GHz. 
Therefore, although the youngest sources (highest gs.4 values) 
have larger ionized gas fractions, there exists a range of source 



compactness for any given age. This causes the lack of correla- 
tion between gs.4 and a l 4 ^. 

5.2. Low frequency radio spectral index and the FIR colour 

As we saw in fig. [7] the low frequency radio spectral index 
a^'g is correlated with the FIR flux ratio /eo/Zioo- Warmer FIR 
colours are therefore associated with more free-free absorption. 
However, we recall that younger sources do not show more free- 
free absorption (Fig. [5] ct\ \ vs. q%.<\). Together, these imply that 
the warmest dust is not found in the youngest sources, as has 
been suggested previously (Vega et al., 2005). If source com- 
pactness has an important role to play in determining the level 
of free-free absorption at low radio frequencies, then the most 
compact are those in which the dust is warmest. 
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5.3. High frequency radio spectral indices 

A small number of ULIRGs in our sample show surprisingly 
steep spectral indices at high frequencies. UGC 5101, Arp 299, 
UGC 8387, UGC 8696, IZw 107, NGC 6286 and NGC 7469 all 
have radio spectral indices between 8.4 and 22.5 GHz steeper 
than a = —I. There appears to be either a deficit in the observed 
flux at 22 GHz or a spectrum that is also very steep at lower 
frequencies (UGC 5101, IZw 107, NGC 6286). In the case of 
IZw 1 07 the radio spectrum has a spectral index of a = -1.4 from 
4.8 to 22.5 GH2Q. In those cases where the radio spectra show a 
steepening towards higher frequencies the spectrum appears to 
show a break around 15 GHz. 

This is contrary to what we expect in sources with a large 
mass of free-free emitting gas. The radio spectrum should flat- 
ten at high frequencies due to the increasing importance of the 
thermal emission (S Y oc y~ 01 ). Although it is conceivable that 
the mass of free-free emitting gas is reduced by the absorption 
of ionizing photons by dust, this effect alone cannot produce an 
observed radio spectrum that steepens towards higher frequen- 
cies. 

A steepening of the radio spectrum towards higher frequen- 
cies is expected in the presence of 'aging' of the relativis- 
tic electron population (Condon, 1992). Because synchrotron 
(and inverse Compton) losses are greater at higher frequencies 
(T sy nc x v~ 5 ) higher energy electrons emit synchrotron pho- 
tons for shorter times than those of lower energies. For an in- 
stantaneous injection of relativistic electrons the radio spectrum 
therefore steepens with increasing frequency. The 'break' fre- 
quency which parameterizes this steepening moves to lower fre- 
quencies at later times after the injection. Spectral aging effects 
were considered by Colbert, Wilson & Bland-Hawthorn (1994) 
to explain the steep radio spectral index in the non-nuclear re- 
gions of NGC 6240. 

In ULIRGs the electron lifetimes are very short because in 
a high radiation energy density environment inverse Compton 
losses would remove all high energy electrons before they could 
emit synchrotron photons. The observation of synchrotron emis- 
sion in such sources therefore shows that they must have strong 
B-fields and short electron lifetimes due to synchrotron losses. 
Electron lifetimes in the GHz region are shorter not only than a 
starburst timescale (10 7 - 10 8 yr) but also than the lifetime of the 
most massive stars (~ 5 x 10 5 yr). Therefore the production of 
relativistic electrons cannot be considered instantaneous, but is 
rather close to a continuous injection approximation. 

Under the assumption of continuous injection of relativis- 
tic electrons, the radio spectrum can be considered as the su- 
perposition of several instantaneously injected populations, each 
with a different age, and therefore each with a different 'break' 
frequency. The spectrum which results from the summation of 
all these components is straight. If the asymptotic change in 
spectral index for a single component of spectral index, a, is 
Aa = -0.5 (such as in the "dynamical halo" model of Lerche 
& schlickeiser, 1981) then the resulting spectrum has a spectral 
index a a g e d = a — 0.5. The radio spectra of all ULIRGs should 
therefore be steepened by the effects of spectral aging, but we 
do not expect to detect a down-turn in the spectra towards high 
frequencies if continuous injection is a good approximation. We 
note however, that because the spectral break for an aging, in- 
stantaneously injected electron population spends little time in 
the ~ 10 GHz range (where we observe the putative breaks) the 



2 Note, however, that a steeper power-law fit to the high frequency 
radio data would not fit the 1 .4 GHz and 408 MHz data points 



form of the time-averaged spectrum is rather sensitive to stochas- 
tic effects in the recent past. Such effects will be important where 
the supernova rate shows variations on timescales of order the 
electron lifetime. 

Qualitatively, such a stochastic effect could be caused by a 
radio hypernova, that might have an order of magnitude greater 
flux than a typical supernova. However, it seems very unlikely 
that a single such source could effect the integrated radio spec- 
trum of an entire ULIRG. The compact sources mapped in 
Arp 220 by Lonsdale et al. (2006) account for only 6% of the 
total flux at 18 cm 

Lisenfeld et al. (2004) considered various explanations for 
the high frequency turn-over seen in the synchrotron spectrum 
of the starburst dwarf galaxy NGC 1569, once the contribution 
of thermal emission had been subtracted. They concluded that 
a sharp turn-over could only be achieved by either a rapid tem- 
poral variation of the star formation rate or through the pref- 
erential escape of low-energy electrons from the galaxy's disk. 
In NGC 1569 a starburst that started abruptly 10 7 yr ago was 
found to reproduce the observed break in the synchrotron spec- 
trum. The escape of low energy electrons from the galaxy can 
produce a break if convection can transport the electrons out of 
the disk during their radiative lifetime. As the lower energy elec- 
trons suffer less synchrotron losses, and therefore emit for longer 
than high-energy electrons, they can in principle be lost prefer- 
entially and cause a break in the spectrum. Lisenfeld et al. (2004) 
found that a convective wind velocity of 150 kms _I can repro- 
duce the spectral break in NGC 1569. However, the lifetimes of 
electrons against radiative losses in ULIRGs should be consider- 
ably shorter than in NGC 1569. Under the assumption of energy 
equipartition between the radiation and magnetic fields, Condon 
et al. (1991a) estimate electron lifetimes of ~ 10 4 yr for elec- 
trons radiating at 8.4 GHz in the same sample of ULIRGs from 
which ours were selected, while that in NGC 1569 was estimated 
to be 5 x 10 6 yr by Lisenfeld et al. Because the supernova rate 
cannot change significantly over timescales of 10 4 yijjand a con- 
vective wind could not remove electrons from the sources over 
similar timescales, none of these mechanisms seem a plausible 
explanation for a high frequency spectral break in the present 
sample of ULIRGs. 

Although there is no doubt that some sources have steep 
spectral indices from 1 to 20 GHz, in the majority of cases where 
a turn-over is seen it is due only to the 22.5 GHz data point 
(Arp 299, UGC 8696, IRAS 15250+3609). We would require 
data at a higher frequency (e.g. 43 GHz) for these sources in or- 
der to be sure of the spectral bend and warrant a more in-depth 
discussion. 



6. Conclusions 

We have presented new high frequency radio data for a sam- 
ple of ULIRGs, that, together with re-reduced archival data and 
fluxes taken from the literature forms a sample of 3 1 ULIRGS 
with well-sampled radio spectra. All but 1 source have a mea- 
sured flux at 22.5 GHz. Every effort has been made to ensure that 
the fluxes are reliable measures of the integrated radio emission. 



3 Type lb and type II supernovae are thought to be the source of 
the synchrotron emission in ULIRGs, and occur by the core-collapse 
of stars more massive than 8 M Q . The difference in lifetime between a 
star of 8 M Q and a star of 100 M Q is ~ 5 x 10 7 yr. Therefore, for a 
normal initial mass function, even if the stars formed instantaneously, 
the supernova rate cannot show large variations on timescales as short 
as 10 4 yr. 
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This means that the possibility of missing flux in high-resolution 
interferometer data and of confusing sources in low resolution 
data has been considered. We find the following, 

- Few sources have straight power-law slopes. Although the 
individual spectra show a variety of spectral bends the trend 
of the mean spectrum is a steepening from 1.4 GHz to 
22.5 GHz. Above and below these frequencies data is too 
sparse to define a meaningful average. 

- The low frequency radio spectral index between 1.4 and 
4.8 GHz is correlated with the FIR-radio flux density ratio 
q calculated at 1 .4 GHz. This correlation is due to free-free 
absorption of low frequency synchrotron photons by ionized 
gas. 

- The high frequency radio spectral index between 8.4 and 
22.5 GHz is correlated with the FIR-radio flux density ratio 
q, calculated at 8.4 GHz. This correlation is due to a higher 
fractional contribution from free-free emission in sources 
with higher values of 4. These are the younger sources. 

- Given the above correlations the lack of a correlation be- 
tween the radio spectral index between 1.4 and 4.8 GHz and 
q%A is evidence that the compactness of the emitting regions 
plays an important role in defining the slope of the low fre- 
quency radio spectrum. 

- The fact that the low frequency radio spectral index a\ \ is 
correlated with the FIR flux ratio /eo A/100 (while a^s K 
not) suggests that the ULIRGs in which the emitting regions 
are most compact are those with the warmest IRAS foal /100 
colours. 

If the FIR-to-radio flux ratio measured at 8.4 GHz, q&A, is 
a measure of the age of the star formation in ULIRGs (larger 
values being 'younger') we find that the youngest objects have a 
larger contribution to their radio fluxes from free-free emission, 
but do not have the most free-free absorption. The sources with 
most free-free absorption are instead those in which the emitting 
regions are more compact. Although these have the warmest FIR 
colours, they are not necessarily the youngest sources. 

Observations at radio frequencies below 1 .4 GHz will be a 
very effective way of quantifying the effects of free-free absorp- 
tion in these sources. 

In paper 2 these radio spectra, together with data from the 
radio to the near-infrared, will be used to constrain new models 
that reproduce the spectral energy distribution of compact star- 
bursts with a possible AGN contribution. 
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Appendix A: Individual objects 

Here we briefly describe aspects of the radio emission relevant 
to the derivation of integrated fluxes. 
IC 1623 

The 22.5 GHz emission shows an unresolved peak at the posi- 
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tion of the eastern nucleus with a clear extension towards the 
west. The distribution of emission is very similar to that of the 
molecular gas (Iono et al. 2004, Yun et al. 1994). 

IIIZw 35 

The radio emission originates from the more northerly of the 2 
optical nuclei. 

IRAS 03359+1523 

1 .'5 to the south of this object there is a 15.4 mJy radio source. 
IRAS 03359+1523 itself emits only 19.4 mJy at 1.4 GHz. Rather 
than reducing the IRAS fluxes by nearly a factor of two we 
choose to exclude this object from the analysis of paper 2. 

UGC 4881 

We detect two sources in this object at the same locations as 
those revealed in the 8.4 GHz observations of Condon et al. 
(1991a). The more easterly component is the brighter and ac- 
counts for ~ 70% of the total flux detected, which is similar to 
the flux ratio at 8. 4 GHz. 

UGC 5101 

A 12.3 mJy radio source lies 1' to the north-east. IRAS fluxes 
reduced by 8%. 

IRAS 10173+0828 

The sub-arcsecond 15 and 22.5 GHz observations of Smith, 
Lonsdale & Lonsdale (1998) show an unresolved nuclear source. 

UGC 6436 

The NVSS shows that a companion galaxy, IC 2810b, located 
1 .'2 to the south-east has a flux of 7.7 mJy compared to the 

19.4 mJy of UGC 6436. This would be confused at IRAS reso- 
lutions and so the IRAS fluxes have been reduced by 26%. No 
other observations were at so low a resolution as to make con- 
tamination from this source a problem. 

Arp 299 

Actually a close interacting pair consisting of NGC 3690 and 
IC 694. 

IRAS 12112+0305 

A source 2 '4 to the south-east has a 1.4 GHz flux of 5.1 mJy 
compared to the 23.3 mJy emitted by IRAS 121 12+0305. IRAS 
fluxes reduced by 15%. 

UGC 8387 

Our 22.5 GHz map shows an unresolved source but the higher 
resolution 8.4 GHz radio continuum map of Condon et al. 
(1991a) shows an elongated structure with an extent of ~ 4". 
Clemens & Alexander (2004) find that free-free absorption flat- 
tens the radio spectral index towards the centre of the source 
showing there to be a large region of dense ionized gas. 

The 22.5 GHz flux for this source may be incorrect. The cal- 
ibrator flux for this source was only about half that reported in 
the VLA calibrator database. 

IZw 107 

Our 8.4 GHz map shows that the extended emission seen at 

22.5 GHz is due to the presence of 2 sources separated by 8". 

Arp 220 

Smith et al. ( 1 998) resolved the compact radio emission into sep- 
arate supernova remnants and subsequent, high sensitivity VLBI 
monitoring of both nuclei has now succeeded in obtaining a di- 
rect estimate of the supernova rate (Lonsdale et al, 2006) of 
4 ± 2 yr _1 . The implied star formation rate is sufficient to supply 
the bolometric luminosity of the system. 

NGC 6286 

The radio structure at 22 and 8 GHz is elongated in the same 



sense as the optical, edge-on disk. A weak source (NGC 6285) is 
present 1 '.5 to the north-west and this is detected in our 8.4 GHz 
data but is undetected in our 22.5 GHz map. The IRAS fluxes 
were reduced by 6% to remove the probable contribution of 
NGC 6285. 

Mrk 331 

See Prouton et al. (2004). 
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Table .1. Radio fluxes. Table .1. cont.d. 



Source 


Frequency 


Flux 


Ref. 


Source 


Frequency 


Flux 


Ref. 




(GHz) 


(mJy) 






GHz 


mJy 




NGC 34 


1.4 


66.9 ±2.5 


d 


UGC 5101 


0.151 


750 ± 60 


f 




4.8 


28.4 ±0.3 


c 




1.4 


170.1 ±5.8 


d 




8.4 


15.2 ± 0.8 


c 




4.8 


73 ± 15 


f 




22.5 


7.41 ±0.14 


b 




8.4 


52.6 


c 


IC 1623 


0.365 


906 ± 83 


V 




15 


30 ± 10 


f 




1.4 


249 ± 9.8 


d 




22.5 


18.0 ± 1 


a 




4.8 


96 ± 12 


i 




85.5 


6 


r 




8.4 


54.6 ± 1 


a 


IR 10173+0828 


1.4 


9.9 ± 0.9 


d 




22.5 


21.7 ±2 


a 




4.8 


5.88 ± 0.2 


c 


CGCG436-30 


1.4 


49.8 ± 1.5 


d 




8.4 


5.4 


c 




4.8 


21.5 


g 




15 


3.1 ±0.3 


c 




8.4 


12.7 ±0.6 


c 




22.5 


2.9 ± 0.4 


a 




22.5 


8.73 ±0.31 


b 


IR 10565+2448 


1.4 


57.0 ± 2.1 


d 


IR 01364-1042 


1.4 


15.8 ± 0.7 


d 




4.8 


22.21 ±0.13 


P 




4.8 


11.8 







8.4 


14.1 


c 




8.4 


8.2 ± 0.4 


c 


Arp 148 


1.4 


36.4 ± 1.2 


d 




22.5 


3.97 ±0.19 


b 




4.8 


14.0 ±0.5 


c 


IIIZw 35 


1.4 


40.6 ± 1.3 


d 




8.4 


11.76 ± 1 


a 




4.8 


26.3 


g 




15 


8.0 ± 1 


c 




8.4 


19.7 


c 




22.5 


6.1 ± 0.3 


a 




15 


10.0 ± 1 


c 


UGC 6436 


1.4 


19.4 ±0.7 


d 




22.5 


9.7 ± 0.25 


b 




8.4 


5.2 ±0.3 


a 


UGC 2369 


1.4 


49 ± 1.5 


d 




22.5 


2.8 ±0.5 


a 




8.4 


13.3 


c 


Arp 299 


0.365 


1550 ±50 


V 




22.5 


5.3 ±0.3 


a 




1.4 


686.3 ± 25.4 


d 


IR 03359+1523 4 " 


1.4 


19.4 ±0.7 


d 




4.8 


403 ± 45 


h 




4.8 


11.7 


g 




8.4 


243 ± 8 


c 




8.4 


11.0 


c 




15 


160 ±3 


c 




22.5 


4.5 ±0.3 


a 




22.5 


74 ± 2 


c 


NGC 1614 


1.4 


137.1 ±4.9 


d 


IR 12112+0305 


1.4 


23.3 ±0.8 


d 




4.8 


63 ± 11 


j 




4.8 


14.7 ±0.2 


c 




8.4 


40.7 ± 0.3 


c 




8.4 


10.0 


c 




22.5 


21.0 ± 1 


a 




15 


5.7 ±0.2 


c 


IR 05189-2524 


1.4 


28.8 ± 1.0 


d 




22.5 


5.7 ±0.3 


a 




4.8 


17.0 ±0.9 


n 


UGC 08058* 


0.365 


551 ± 39 


V 




8.4 


11.4 


c 




1.4 


274 ± 3 


c 




15 


7.8 ±0.2 


c 




4.8 


265 ± 1 


c 


NGC 2623 


1.4 


95.7 ±2.9 


d 




8.4 


189 ± 1 


c 




2.4 


83.0 ±4 


1 




15 


146 ± 1 


c 




4.8 


59.0 ± 10 


h 




22.5 


137 ± 2 


c 




8.4 


35.5 


c 


UGC 8387 


1.4 


104.4 ± 3.2 


d 




22.5 


18.3 ±0.2 


a 




4.8 


46 


m 


IR 08572+3915 


1.4 


4.3 ±0.4 


d 




8.4 


34.9 


c 




4.8 


4.4 ±0.2 


n 




15 


12 ± 2 


c 




8.4 


4.1 ±0.2 


c 




22.5 


7.5 ± 2 


a 




22.5 


3.18 ±0.30 


b 


NGC 5256 


0.365 


353 ± 43 


V 


UGC 4881 


1.4 


36.8 ± 1.2 


d 




1.4 


126.3 ± 4.5 


d 




8.4 


8.8 


c 




-1 o 

4.8 


43.3 ± 3.1 


n 




22.5 


3.5 ±0.5 


a 




8.4 


24.1 ± 1 


c 


+ Only one source of a close pair emits in the radio. * Source possibly 




15 


21.0 ± 1 


c 


variable in radio. 1.4, 4.8, 


8.4, 15 and 22.5 GHz fluxes all made on the 




20.0 


14.9 ± 1.5 


n 


same date. References, (a) this work, (b) Prouton et al. (2004). (c) 


UGC 8696 


0.151 


590 ± 60 


f 


reduced from VLA archive. See table[3]for details, (d) NVSS, Condon 




0.365 


335 ± 37 


V 


et al. (1998). (e) Condon et al. (1991a). (f) Sopp & Alexander (1991). 




1.4 


144.7 ± 5.1 


d 


(g) Sopp & Alexander (1992). (h) Gregory & Condon (1991). (i) 




2.7 


96 ± 6 


q 


Griffith et al. (1994). (j) Griffith et al.(1995). (k) Condon et al.(1983). 




4. 8 


71 + 15 


f 


(1) Dressel & Condon (1978). (m) Condon, Frayer & Broderick 




8.4 


43.5 


c 


(1991b) . (n) Rush, Malkan & Edelson (1996). (o) Baan & Klockner 




15 


29.5 ± 1 


c 


(2006). (p) Crawford et al. (1996). (q) Sramek & Tovmassian (1976). 




22.5 


14.3 ±0.3 


a 


(r) Imanishi, Nakanishi & Kohno (2006). (s) Zhao et al. (1996). (t) 




85.6 


6 


r 


Rodriguez-Rico et al. (2005). (u) Anantharamaiah et al. (2000). (v) 


IR 14348-1447 


1.4 


35.9 ± 1.2 


d 


Douglas et al. (1996) (w) Ficarra, Grueff & Tomassetti (1985) 




8.4 


10.8 ± 1 


a 












22.5 


4.5 ± 0.5 


a 
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Table .1. cont.d. 



Source 


Frequency 


Flux 


Ref. 




GHz 


mjy 




IZw 107 


0.408 


100 ± 20 


w 




1.4 


50.3 ± 1.6 


d 




4.8 


32.0 ± 0.6 


h 




8.4 


14.3 ± 1 


a 




22.5 


3.6 ±0.3 


a 


IR 15250+3609 


0.151 


< 100 


f 




1.4 


14.5 ±0.6 


d 




4.8 


13 ±3 


f 




8.4 


10.5 


e 




15 


8±3 


f 




22.5 


4.2 ± 0.3 


a 


Arp220 


0.151 


260 ± 30 


f 




0.365 


435 ± 26 


V 




1.4 


326.3 ± 9.8 


d 




2.4 


312 ± 16 


1 




2.7 


260 ± 13 


k 




4.8 


210 ± 20 


f 




8.4 


148.0 


e 




15 


110 ± 1 


f 




22.5 


90 ±6 


s 




43 


44 ±4 


t 




97.2 


61 ± 10 


u 


NGC 6286 


0.365 


431 ±37 


V 




1.4 


156.7 ±5.6 


d 




4.8 


53.0 ±8 


h 




8.4 


31.5 ± 1.6 


a 




22.5 


11.2 + 1.5 


a 


NGC 7469 


0.365 


397 ± 24 


V 




1.4 


180.5 ± 5.4 


d 




4.8 


70 ±8 


g 




8.4 


50.2 ± 0.6 


c 




22.5 


17.50 ± 0.5 


b 


IC 5298 


1.4 


34.7 ± 1.4 


d 




4.8 


13.9 


g 




8.4 


8.2 ± 0.4 


e 




22.5 


3.86 ±0.18 


b 


Mrk331 


1.4 


70.7 ± 2.2 


d 




4.8 


31.1 


g 




8.4 


22.3 ± 1 


a 




22.5 


10.0 ± 0.3 


b 
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Fig. 3. Contour maps for the archival VLA data. In brackets after the name of the source is the name of the radio frequency band; 
L (1.4 GHz), C (4.8 GHz), X (8.4 GHz), U (15 GHz), K (22.5 GHz). Contour levels are 2 " c ' /2 where x = 0, 1, 2, . . .The dotted 
contour has the same absolute value as the first contour but is negative. The value of n for each contour map is included in table [3] 
The synthesized beam is shown in the bottom left of each panel. 
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Fig. 4. The radio spectra of the sample galaxies from 100 MHz to 100 GHz. The lower left panel shows the mean spectrum of all 
the sources scaled so that the 1.4 GHz flux is 0.1 Jy. 
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Fig. 5. Relationship between radio spectral index at various frequencies and the logarithmic FIR-radio flux density ratio q, defined 
as q = log(FIR/3.75 x 10 12 Hz)/S, where S is the radio flux in units ofWm^Hz" 1 . The top panels show the radio spectral index as 
a function of q measured at 1.4 GHz and the bottom panels show the radio spectral indices as a function of q measured at 8.4 GHz. 
The left-most panels show the radio spectral index over the whole range from 1 .4 to 22.5 GHz while the remainder show the spectral 
indices over a more restricted range. 



